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Abstract 

We present an exact analytic study on the topological phase interference 
effect in resonant quantum tunneling of the magnetization between degener- 
ate excited levels for biaxial ferromagnets in an arbitrarily directed magnetic 
field. We show that the topological phase interference effect depends on the 
orientation of the field distinctly. The transition from classical to quantum 
behavior is also discussed. 
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During the last decade, there has been great interest in the problem of quantum tunnel- 
ing in nanometer-scale magnets.0 One notable subject is that the topological Wess-Zumino 
(or Berry) phaseB can lead to remarkable spin-parity effects. It was theoretically shown 
that the tunnel splitting is suppressed to zero for half-integer total spins in single-domain 
biaxial ferromagnetic (FM) particles due to the destructive interference of the Berry phase 
between two tunneling paths of opposite windings.il However, the interference is constructive 
for integer spins, and hence the splitting is nonzero. While spin-parity effects are sometimes 
be related to the Kramers degeneracy, they typically go beyond the Kramers theorem in a 
rather unexpected way.SH The auxiliary particle method was proposed to study the spin- 
parity effects in one model of a single large spin subject to the external and anisotropy 
fields.i Recently, the spin-phase interference and quantum oscillation effects were studied 
extensively in FM particles in the presence of a magnetic field, with the field along either 
the hard,B0 easy,i or medium axis.i Similar spin-parity effects were also found in antiferro- 
magnetic (AFM) particles,00 and in the quantum propagation of Bloch walls in quasi-one- 
dimensional ferromagnets0 and antiferromagnets.lii By applying the effective Hamiltonian 
approach, the effect of magnetic field and quantum interference of the magnetization vector 
(or the Neel vector) were studied in FM (or AFM) particles with different symmetry.0 

Experiment on the molecules Fes showed a direct evidence of the role of the topological 
phase in spin dynamics.El Recent theoretical and experimental studies include the spin 
tunneling in a swept magnetic field,0 the thermally activated resonant tunneling based on 
the perturbation theory0 and the exact diagonalizationjli the discrete WKB method and 
a nonperturbation calculationjll the non-adiabatic Landau-Zener model,@ the calculation 
based on exact spin-coordinate correspondence, ill and the effects caused by the higher order 
term and the nuclear spins on the tunnel splitting of Feg.i^ 

It is noted that the previous studies on FM spin-parity effects!' Wi were mostly focused 
on the phase interference between two opposite winding ground-state tunneling paths. More- 
over, the previous work^UH have been confined to the condition that the magnetic field 
be applied along the easy, medium, or hard axis, separately. In this paper we study the 



topological phase interference effects for single- domain biaxial FM particles in an arbitrar- 
ily directed magnetic field. Our study provides a nontrivial generalization of the Kramers 
degeneracy for equivalent double-well system to coherently spin tunneling at ground states 
as well as low-lying excited states for FM system with asymmetric twin harriers caused by 
the arbitrarily directed magnetic field. Integrating out the momentum in the path integral, 
the spin tunneling problem is mapped onto a particle moving problem in one-dimensional 
periodic potential U (0). By applying the periodic instanton method, we obtain exactly the 
splittings between degenerate excited levels of neighboring wells. U (0) is regarded as a one- 
dimensional superlattice with a periodically recurring asymmetric twin barriers. The general 
translation symmetry results in the energy band structure, and the low-lying energy level 
spectrum is obtained by using the Bloch theorem. Our result shows that the excited-level 
splittings depend significantly on the parity of total spins of FM particles, and this spin- 
parity effect depends on the orientation of the field distinctly. The transition from quantum 
to classical behavior is also studied, and the second-order phase transition is shown. 

For a spin tunneling problem, the tunneling rate is determined by the Euclidean transi- 
tion amplitude 



and 5* is the total spins. The first two terms in Eq. (2) define the topological Wess-Zumino 
term which arises from the nonorthogonality of spin coherent states.il The total derivative 
has no effect on the classical equations of motion, but is crucial for the spin-parity effects. 

The system of interest has the biaxial symmetry, with x being the easy axis, y being the 
medium axis, and z being the hard axis. The magnetic field is applied in the ZY plane, at 
an angle in the range of | < 6'// < tt. Now E {6, (p) in Eq. (2) is 




(1) 



where DQ = sin OdOdcj), and the Euclidean action is 




MqV = Ai V = h'-fS, where V is the volume of the particle, 7 is the gyromagnetic ratio. 
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E {0, (f)) = K±_ cos^ e + K\\ sin^ 9 sin^ - Mq/Zj, sin ^ sin - MqH^ cos (3) 

where Hy — HsmOff, — HcosOh, K\\ and K^^ are the longitudinal and the transverse 
anisotropy coefficients satisfying K±_ » K\\ > 0. As K± ^ i^y > 0, the deviations of 9 
about 6*0 are small. Introducing 6* = 6*0 + 0;, la] <^ 1, Eq. (3) reduces to 

E {a,(j)) = K± sin^ 6oa^ + K\\ sin^ 60 (sin — sin 0o)^ , (4) 

where cos^o = sin0o = = , X = h = ^, and i/o = 

Performing the Gaussian integration over a, we obtain the transition amplitude as 



Ke — exp [—iS (1 — cos^o) (4>f — J d4>^W \~ J \^ ~'~ ^ 



(</>) 



(5) 



with m — 2^^v ' ~ -f^ii^sin^ ^0 (sin0 — sin^o)^- Since the configuration space of 

this problem is a circle, calculations are restricted to the first twin barrier at = | and 
We use A to denote the instanton passing through the small barrier at = | with the height 
hUs — K\\Vsm'^ 9q (1 — sin0o)^, and B through the large barrier at = ^ with the height 
UUl = i^ii^ sin^^o (1 + sin0o)^- Correspondingly, there are two kinds of anti-instantons: 
A- and B' . 

The periodic instanton configuration 0p which minimizes the Euclidean action in Eq. (5) 
at an energy E > Q satisfies the equation of motion 

Then we obtain the periodic instanton A solution as 

. l-{isn2(a;ir,/ci) 

corresponding to the transition of (f) from arcsina to tt — arcsina, where a = sin0o + 
■\J K\\V^n^ - sn{ijJiT,ki) is the Jacobian elliptic sine function of modulus ki, where kf = 

{iriyiSg' = ^^<i>o - Jj^Jg^,, 6 = i^^, = ff, = 2^^^^sineo, and g. 
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V(l+a)(l-/9) 



. The classical action or the WKB exponent in the tunnel splitting is obtained 



by integrating the Euclidean action in Eq. (5) with the above periodic instanton solution. 
The result for instanton A is Sa = Wa + ^E/S, where 



Wa = 4:1710; I 



(8) 



K (ki), E (ki), and 11 {ki, ^i) are the complete elliptic integral of the first, second, and third 
kind. 

Now we discuss briefly the calculation of the tunnel splitting. For a particle with mass m 
moving in a smooth double- well-like one-dimensional potential U (x), the instanton approach 
gives the ground-state tunnel splitting aJ^l 

AEo = 2C(-^) exp{-Wo), (9) 



27r 

where Wq is the classical action for the ground-state tunneling, 

C f detM^+g;^) 

\a^t'[-d? + U"{x,i{T))/m]j ^ ^ 

is the ratio of fluctuation determinants, and Xd (t) is the instanton solution. The prime on 
the det indicates that the zero eigenvalue is to omitted, and u is the frequency of harmonic 
amplitude oscillations in the wells about the minima. The general formulas were presented 
to evaluate the preexponential factor in the tunnel splitting or the tunneling ratciliil It was 
found that what is need for this evaluation is the asymptotic (r ±oo) behavior of the 
instanton velocity, 

dxci 



dr 



aexp (=F^t) , as r — * ±oo. (11) 



Then the ground-state tunnel splitting isc 



J24 







a 





1/2 

exp i-Wo) . (12) 

In most physical applications, however, the ground-state tunnel splitting can be best esti- 
mated as 

AEo=PoUj[^^) exp{-Wo), (13) 



where the dimensionless prefactor po is often relevant to experiments. It is noted that 
Eq. (12) is based on quantum tunnehng at the level of ground state, and the temperature 
dependence of the tunneling frequency (i.e., tunneling at excited levels) is not taken into 
account. The instanton technique is suitable only for the evaluation of the tunneling rate 
or the tunnel splitting at the vacuum level, since the usual (vacuum) instantons satisfy 
the vacuum boundary conditions. Recently, different types of pseudoparticle configurations 
(periodic or nonvacuum instantons) are found which satisfy periodic boundary conditions.!^ 
For the same tunneling problem, the WKB answer for the splittings of the nth degenerate 
excited levels or the imaginary parts of the metastable levels idHi@ 

AEn (or Im E„) = '^^-^ exp {-W) , (14) 

vr 

where u (En) = is the energy-dependent frequency, t (_E„) is the period of the real-time 
oscillation in the potential well, 

t (En) = V2^ / , (15) 

where Xi^2 {En) are the turning points for the particle oscillating inside the potential U (x). 
The functional-integral and the WKB method showed that for the potentials parabolic near 
the bottom the result Eq. (14) should be multiplied by *'^"'2"e"n! ^ '^^i?> factor is very 
close to 1 for all n: 1.075 for n = 0, 1.028 for n = 1, 1.017 for n = 2, etc. Stirling's formula 
for n\ shows that this factor trends to 1 as n oo. Therefore, this correction factor, 
however, does not change much in front of the exponentially small action term in Eq. (14). 
It is noted that Eq. (14) can be obtained by using the connection formulas near a linear 
turning point ,@ or by matching the WKB wavefunction in the classically forbidden region to 
the exact harmonic oscillator wavefunction near the classical turning point .0 It was shown 
that the WKB method is equivalent to the instanton method for the ground-state tunnel 
splittingB Liang et al. showed that the WKB method and the periodic instanton method 
give the same result for the excited-state tunnel splitting.il 

For the present problem, we find that As a = ^^^^^ exp {—Wa), where tA^E) = 

5 



^)K(k[) and k[ = ^ — kf. The same method can be apphed to the instanton B passing 
through the large barrier. And the result is Aeb = j^e) ^-^P for < < Us, where 

ts (E) — tA (E), and Wb has the same expression as Eq. (8) but taking as ^2 = jz^- 
For Us < E < Ul, the imaginary part of the metastable level is Im£^ = ^tt^t exp (—2Wb), 

tB\E) \ J 

where iB{E) = ^^K{k'^), 



Wb = 2m^3 (1 + a) 002 



-2 — -E {k2) - -K (/cs) + n (A;2, Cs) 



(16) 



with k'2 = Jl - kl UJ2 = ^,92= ^/S' ^3 = and k^ = ^Z'af^ - Now we discuss the 



low energy limit of the level splitting. With the help of harmonic oscillator approximation 
for energy near the bottom of the potential well, En = (n + | j il, and 



1 fdW\ 12K,V . ^ 

■ sm cos 00) 



'\ m \ o?0^ J ^^^^ V hm 



WAiB),n = WAiB),o - ( n + ^) + (n + ^) In (^^^±11^-—] , (17a) 



Wa{b) can be expanded as 

2/ V 2/ \8V A-Ssin^ocos^ 00/ 
W^A(s),o = 2\/XS sin 6*0 ^cos 0o -F 2 sin 0o arcsin 2 ' (■^'^^) 

where "— " for the instanton A, and "+" for the instanton B. Therefore, the tunnel splittings 
for instantons A and B are 

hAsA{B),n = (-^-L^) V^-S""-^ sin 610 COS 00 (Sx/A^ sin 9o cos^ 0o) ^ exp (^-Wa{b),o) ■ 

(18) 

Equation (17b) shows that the WKB exponent for instanton A is smaller than that for 
instanton B at finite magnetic field because the barrier through which instanton B must 
tunnel is higher than that for instanton A. For the case of ground-state resonance, i.e., 
n — 0, Eq. (18) reduces to 

A£a(s),o = (sin ^0 cos^ 0o) ^ exp (-W^A(i?),o) • (19a) 



It is not particularly illuminating to write out the prefactor pq (see Eq. (13)), although it 
can be seen that it is dimensionless and independent of the volume V of the particle. Note 
that Eq. (13) is the approximate formula for ground-state tunnel splitting with exponential 
accuracy. In most physical applications, the preexponential factor can be best estimated as 
an attempt frequency. The apparent disagreement with Eq. (13) was also found in other 
spin-tunneling problems.^ However, in the case of zero magnetic field Eq. (19a) reduces to 



AsAfl iH = 0) = AsBfl (H = 0) 



4^0 ( ^ ) exp i-Wo) 



(19b) 



where = ^ (^ = 0), and Wo = Wa,o {H = 0) = Wb,o {H = 0) = 2^S. Compared with 
Eq. (13), po = 4 for this case. 

It is noted that TiA.eA{B),n is only the level shift induced by tunneling between degen- 
erate excited states through a single barrier. U {(p) can be regarded as a one-dimensional 
superlattice with the sublattices A and B. The Bloch states for sublattices A and B are 
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and 



(e,0) = ^Ee'^'^"^A (0-0n), 



$5 (e, 0) = 4f E e'^(^"+"ViJ (0 - 0n - a) 



(20a) 



(20b) 



where 0n = Iwu + (po, L = N {a + b) , a = it — 20o, and 6 = tt + 20o- The total wavefunction 
is 



^5 (0) = aA (0 {^,<P) + aB (0 (e,' 



(21) 



Including the Wess-Zumino phase, we derive the secular equation in the tight-binding ap- 
proximation as 



OA (0 

as (0 



0, (22) 
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where fj, — S {1 — cos^o), and the Bloch wave vector ^ = in the first Brillouin zone. Then 
the tunnel sphtting of the nth excited level is 



Aen = 2sJ{AeA,nf + {^eB,nf + 2 {AeA,n) (A£b,„) cos [2t:S (1 - cos^o)]- (23) 

Another approach to obtain Eq. (23) is to calculate the transition amphtude directly. The 
subtle point in evaluating the transition amplitude is how to arrange the instantons and 
anti-instantons appropriately to satisfy the boundary condition. Note that one configuration 
starting from |^o,0o) and ending at |^o,0o) can be an arbitrary permutation of the pairs 
{AB), {AA~), {B~B), and {B~A~). Introducing s, t, p, and q as the numbers of A, B, 
A~ , and B~ in the instantons and anti-instantons pairs, and i, j, k, and / as those of {AB), 
{AA~), (B-B), and {B-A'), we have 

ijf-j^s,i + k^tj + l^p,k + l^q. (24) 

Therefore, s + q — t + p, and only three variables are independent. Now the transition 
amplitude is 

X^^,&-("+^)" V N{s,t,p,q) Wf^^ TY+P^f^^ Ty+i 

(25) 



^g-iS(l-coseo)(7r-2(/.o)(s-p)g-iS'(l-cos 6»o)(7r+2(/.o)(t-g) 



where N {s,t,p,q) represents the number of different configurations for a given set of 
{s, t,p, q}. It can be calculated as 

i=max{0,t—q} 



which has a simple result 



Then we have 

Ke = y^e-("+^)^cosh y {AeA,nf + {AeB,nf + 2 {AeA,n) ( A£b,„) cos [2nS (1 - cos ^0)]^ 

(28) 



and we can read off the splitting of the nth excited level shown in Eq. (23) directly. 

Finally we discuss the phase transition from classical to quantum behavior. It was found 
that for a particle moving in a double- well potential U (x), the behavior of the energy- 
dependent period of oscillations P {E) in the Euclidean potential —U {x) determines the 
order of the quantum-classical transition.ii If P (E) monotonically increases with the am- 
plitude of oscillations, i.e., with decreasing energy E, the transition is of second order. The 



corresponding crossover temperature is Tg = ujq = y ^ \U" {xsad)\ is the frequency of 
small oscillations near the bottom oi —U (x), where Xsad corresponds to the top (the saddle 
point) of the barrier. If, however, the dependence of P (E) is non-monotonic, the first order 
crossover takes place. For this case, the period of the instanton A is Pa{E) = ^K{ki). 
The monotonically decreasing behavior of Pa (E) in the domain < E < Us is shown in 
Fig. 1 for several values of sin0o5 which shows that the second-order phase transition takes 
place, and the crossover temperature is 



The period of the instanton B is Pb (E) = Pa (E) for < E < Us, while Pb (E) = (fcg) 
for Us < E < Ul- Fig. 2 shows the monotonically decreasing behavior of Pb (E) for the 
whole domain < E < Ul- Again one finds a second-order transition from the thermal to 
quantum regime, and the crossover temperature is 



At zero magnetic field, AeA,n = ^^B,n, the tunnel splitting is suppressed to zero for the 
half-integer total spins by the destructive interfering Wess-Zumino-Berry phases, which is in 
good agreement with the Kramers theorem. The presence of a magnetic field perpendicular 
to the plane of rotation of magnetization yields an additional contribution to the Berry 
phase, resulting constructive and destructive interferences alternatively for both integer and 
half-integer spins. Tunneling is thus periodically suppressed. At finite magnetic field, the 
low-lying tunneling level spectrum depends on the parity of total spins. If S is an integer. 





(29a) 




(29b) 
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the tunnel splitting is 



A£„ = 2^(A£A,n)' + {^eB,nf + 2 (As^.n) (A£b,„) COS (27r,S COS ^o) ■ 



(30a) 



While if 5" is a half-integer, the sphtting is 



A£„ = 2^(A£A,n)' + (A£b,„)' - 2 (As^.n) (A£b,„) COS (27r5cos^o)- 



(30b) 



Equations (30a) and (30b) clearly show oscillations of the tunnel splitting as a function of the 
transverse magnetic field, together with the parity effect. Note that the magnetic field along 
the medium axis does not produce any oscillations. As ^£A,n > Ae^.n at finite magnetic 
field, the tunnel splitting will not be suppressed to zero even if the total spin is a half-integer. 
The similar oscillation of tunnel sphtting with the magnetic field was also found in Ref. 29, 
while the system considered in Ref. 29 is a single-domain antiferromagnetic particle with 
a field along the hard axis. The most interesting observation is that only the z component 
of the magnetic field (i.e., along the hard axis) can lead to the oscillation of the tunnel 
splitting. As shown in Fig. 3, the splitting depends on the magnitude of the field H and its 
angle 6h with the hard axis, and even a small misalignment of the field with the z axis can 
completely destroy the oscillation effect. For small Oh the tunnel splitting oscillates with 
the field, whereas no oscillation is shown up for large 6h- In the latter case, a much stronger 
increase of tunnel splitting with the field is shown. This strong dependence on the orientation 
of the field can be observed for ground-state resonance as well as excited-state resonance. 
As a result, we conclude that the topological phase interference or spin-parity effects depend 
on the orientation of the external magnetic field distinctly. The distinct angular dependence, 
together with the oscillation of the tunnel splittings with the field, may provide an independent 
experimental test for the spin phase interference effects in FM particles. This "Aharonov- 
Bohm" type of oscillation in magnetic system is analogous to the oscillations as a function of 
external flux in a SQUID ring. Due to the topological nature of the Berry phase, these spin- 
parity effects are independent of details such as the magnitude of total spins and the shape 
of the soliton. The transition from classical to quantum behavior is also studied for the small 
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and large tunneling barrier, respectively. By calculating the periods of instanton A and B 
analytically, we find the monotonically decreasing behavior of the periods with increasing 
energy, which yields a second-order phase transition. Our results should be useful for a 
quantitative understanding on the spin-parity effects and the quantum-classical crossover 
for FM particles in an arbitrarily directed magnetic field. It is noted that the theoretical 
results presented here are based on the instanton method, which is semiclassical in nature, 
i.e., valid for large spins and in the continuum limit. The theoretical calculations performed 
in this paper can be extended to the single-domain antiferromagnetic particles where the 
relevant quantity is the excess spin due to the small noncompensation of two sublattices. 
Work along this hne is still in progress. We hope that the theoretical results obtained in the 
present work will stimulate more experiments whose aim is observing the topological phase 
interference effects and the quantum-classical transition in resonant quantum tunneling of 
magnetization in nanometer-scale single-domain ferromagnets. 
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Figure Captions: 

Fig. 1 The relative period Pa (E) /Pa{E — Us) of the periodic instanton A as a function 
of energy E/K\\V sin^ 9q in the domain < < Us- 

Fig. 2 The relative period Pb {E) /Pb{E — Us) of the periodic instanton S as a function 
of energy E/K\\y 9o in the domain < £■ < Ul- 

Fig. 3 The relative tunnel splitting of the ground level (n = 0) Aso/Asb q {H = 0) as a 
function of H/Hq for Oh = 0°, 1°, 3°, 5°, and 90°. Here A = K\\/K± = 0.1. 
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